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Abstract 

We propose a self-consistent approximate solution of the disordered Kondo-lattice model (KLM) 
to get the interconnected electronic and magnetic properties of 'local-moment' systems like diluted 
ferromagnetic semiconductors. Aiming at [Ai-xMx) compounds, where magnetic (M) and non- 
magnetic (A) atoms distributed randomly over a crystal lattice, we present a theory which treats the 
subsystems of itinerant charge carriers and localized magnetic moments in a homologous manner. 
The coupling between the localized moments due to the itinerant electrons (holes) is treated by a 
modified RKKY-theory which maps the KLM onto an effective Heisenberg model. The exchange 
integrals turn out to be functionals of the electronic selfenergy guaranteeing selfconsistency of our 
theory. The disordered electronic and magnetic moment systems are both treated by CPA-type 
methods. 

We discuss in detail the dependencies of the key-terms such as the long range and oscillating 
effectice exchange integrals, 'the local-moment' magnetization, the electron spin polarization, the 
Curie temperature as well as the electronic and magnonic quasiparticle densities of states on the 
concentration x of magnetic ions, the carrier concentration n, the exchange coupling J, and the 
temperature. The shape and the effective range of the exchange integrals turn out to be strongly 
x-dependent. The disorder causes anomalies in the spin spectrum especially in the low-dilution 
regime, which are not observed in the mean field approximation. 

PACS numbers: 75.50.Pp, 71.23.-k, 75.30.Hx, 85.75.-d 
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I. INTRODUCTION 



There are many analytical models such as the Hubbard and the Anderson models which 
are very useful for the description of real correlated electron systems. The Kondo-lattice 
Model (KLM) is yet another one. The KLM describes an interplay of itinerant electrons in a 
partially filled energy band with magnetic moments localized at certain lattice sites 

naa. 

The characteristic model properties result from an interband exchange interaction between 
two well-defined subsystems: itinerant electrons and localized spins. This model has a 
long history, particularly concerning with itinerant electron magnetism ^, heavy fermion 
systems jsl, perovskite manganese oxides and so on. 

Problems connected with the substitutional disorder have recently become more and more 
important in different fields of material science, e.g., the diluted magnetic semiconductors 
(DMS) j^, the transition metal dielectrics [3], etc. Extensive theoretical work on the problem 
of disorder has been done for random Heisenberg spin systems 8,9,10,[ll,12,13,14, Is], 16 1. 

Here, we are mainly interested on how the disorder influences the characteristic properties 
of local-moment systems such as the DMS, where magnetic (M) and non-magnetic {A) 
atoms are distributed randomly over a crystal lattice {Ai_xMx) with a given concentration 
of rnagnetic atoms x. In order to answer this question different approaches were proposed 8|, 
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16 



20 



25 



261]. More or less realistic electronic structure calculations based on density 



functional theory (DFT) or numerical methods like classical Quantum Monte Carlo were 
used for standard models like the Heisenberg spin system . However, the disadvantage of 
the realistic DFT-calculations come from their strong material dependence 25||. Therefore, 
they do not explain the basic physics of disordered local-moment systems in a simple way. 
But, it is known that a disordered KLM Hamiltonian can be mapped onto an random 
Heisenberg model(shown in the text), so, in this respect better results can be obtained from 
its model study. 

A special challenge when treating the random Kondo-lattice model arises with the fact 
that both the electron and the spin subsystem have to be considered simultaneously and 



on the same level. Most of the KLM investigations are focused on the electronic 28,l34| or 



magnetic 
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33| subsystem only. A special goal of our study is the homologous 



treatment of the electronic and magnetic properties of the random KLM, which mutually 
effect each other and, therefore, should be determined self-consistently. 
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Nevertheless, we are forced to apply different methods to study the influence of the dis- 
order and dilution of the magnetic moments subsystem on the properties of these two afore- 
mentioned subsystems. In this paper for the itinerant electron system a proper alloy analogy 
with the respective coherent potential approximation (CPA) 24| is used. For the random 
spin system, for which the situation is not so clear, an equivalent ansatz must be found. 
As it was done successfully for the periodic KLM []]] {'modified' RKKY (MRKKY)) one 
can map the KLM-interband exchange on an effective and random Heisenberg model. The 
resulting effective exchange integrals between the localized spins will be long range and com- 
plicated functionals of the electronic self-energy. The conventional RKKY, resulting from 
second order perturbation theory is insufficient even with a phenomenological damping fac- 



tor 



3, 



331]. Higher order conduction electron self-energy effects, being taken into account by 



'modified' RKKY but neglected by conventional RKKY, provide the self-consistency of the 
full KLM. They drastically influence the magnetic properties such as the Curie temperature. 

A similar strategy of treating the disordered Kondo-lattice model has recently been used 
in ref. 35|], namely the combination of partial solutions for the electronic and magnetic part 
by a modified RKKY theory. However, the authors use different methods in order to solve 
the partial problems. For the electronic part, they apply a procedure which gives correct 
results only for very low band occupation. For the spin part, on the other hand, the disorder 
has been solved by a mean-field ansatz. In this paper we present alternative treatments of 
these partial problems. In particular, we develop a better technique to deal with the spin 
disorder. The band electrons are discussed within CPA starting from an atomic limit alloy 
analogy without any restriction to the electron band occupation. Several anomalies in the 
spin spectrum are found, which are not reproducible by the method of ref. 35 1. 

The paper is organized as follows: In Sec. 2 we describe the methods and approximations 
that we have chosen to study electron and spin dynamics in the disordered KLM. Numerical 
results concerning Curie temperatures, the magnon densities of states and the effective 
exchange integrals are presented in Sec. 3 for wide range of the model parameters like 
exchange coupling J, concentration of magnetic atoms x, electron band occupation n, and 
temperature T. Section 4 is assigned for a summary and an outlook. 
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II. THEORETICAL MODEL 



The correlated KLM-Hamiltonian can be written in second quantized form as the sum of 
a kinetic energy, an exchange interaction, and a Hubbard-type Coulomb interaction 

= ''^^tijaf^aja- + — '^^nifyTii^^cr — {^o-'S'j^fli^Ojo- + •S'f a^„^ajo- } , (1) 

where; 

z^ = 6^^-6^^,S^ = S! + tz^Sl (2) 

and af^ (cLia) is the creation (annihilation) operator for the Wannier electron with the spin 
a (cr =1, 1) at the site Ri, J is an exchange coupling which we assume to be positive in the 
following theory, tij is an electron hopping integral between lattice sites, U is the Coulomb 
interaction, and is the total number of crystal sites. The inclusion of the Hubbard-term 
to the standard-KLM has only the reason to push doubly occupied levels above the Fermi 
level. This will be further commented on in section 2.1 where we discuss the electronic 
subsystem. 

In order to introduce the disorder into KLM of the Ai_xMx type, we introduce projection 
operators; 

Xf 

Using these definitions we can change the notations in the Hamiltonian ([1]) for the disorder 
in hopping and exchange coupling terms, respectively; 

J. xkk vk -yk 

7 , jk yk 

l^k&A,M 1 

where hopping parameters t"^^ = t'^^^^ = t^^^ = t^^^ = t = W/6 are the same for the 
different kind of atoms and are equivalent to a half band width W for the case of a simple 
cubic lattice, and exchange couplings are J*^ = J, J"^ = for magnetic and nonmagnetic 
atoms, respectively. 



1, - if site i is A 

(3) 

0, - otherwise, 



1, -if site i is M 

(4) 

0, - otherwise. 
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FIG. 1: General procedure of the self-consistent mapping of the disordered Kondo-lattice model 
onto disordered Heisenberg model 

After configurational averaging these operators are simply expressed by the concentration 
X of magnetic atoms M. 

IXf) = l-x, 

\ Wc g 

(xn. = - 

In order to study conduction-electron properties, we use the configurationally averaged 
single-electron Green function; 

G.,. {E) = ((a..|a+»^. (7) 

Where the symbol (...) denotes the configurational ensemble average. 

The relation between the band occupation n and the chemical potential fi is as follows: 

n = n^ + n^ = --}^ (8) 

with j3 = 1/kT is the inverse temperature. 

The main idea of this paper is to map the disordered KLM on an effective random 
isotropic spin Heisenberg Hamiltonian; 

where Jij is an effective exchange interaction between localized magnetic moments. We 
assume that we can calculate the effective exchange interaction, J^j, using the MRKKY 
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method 

m 

j2 poo rl P 1 r 1 

= L i«E^7v ? N ''f-.> C"] • 

which in the weak-couphng hmit agrees with conventional RKKY, where is the Fourier- 
transform of Jij. Then we get the full self-consistent loop for a finite temperature calculation 
(see Fig. [1]). In order to study magnetic properties of the disordered KLM we calculate the 
configurationally averaged magnon Green function: 



D,iE) = {{S^\SZ,))^. (11) 

Using Callen equation 



3l| . it is then possible to calculate the magnetization (5"^): 

{S - 0) (1 + 0)'^+' + {S+1 ^ ^25+1 
1 + 0)^^+^ - 025+1 



where 

The procedure presented in Fig. [1] is very general. It is clear that we have to use different 
approximations in order to solve the disordered electronic and magnetic parts. Now, we 
discuss these approximations in more detail. 

A. Electron Subsystem: Zero-bandwidth limit of the correlated KLM 

The total Hamiltonian of the correlated KLM model is given by ([!]). 

n 

The zero-bandwidth limit [Iq] is defined by 

- 0. (14) 
In this approximation the excitation spectrum consists of the following four poles |l^ 

e, = -lJS, e2 = lJ{S + l), ^^^^ 
e^ = U+\JS, t^ = U -\J{S + 1). 

It means that the single electron spectral density A„ {E) must be a four-pole function 

4 

(E) = J2 {E - e„) . (16) 

m=l 



The temperature- and concentration-dependent coefficients ama have the physical meaning 
of spectral weights for the corresponding excitation energies. The expressions for these 
weight-factors are 



where A^- = (^S°'at(ja„) + {S^ria) is a mixed spin-electron correlation function. 
The term A^ can be expressed by the single-electron Green function [17j : 



(17) 



1 1 rIF' 
Aa = - — — V/ -^7w^, (E~tAlmGr (E). (18) 

k 

A propagating a electron will meet at a certain lattice site Ri the atomic level ei with 
probability ai^, the level €2 with probability a2cr and so on, if there is no correlation between 
sites. This leads to the four- component alloy. Eq. f|T5|) makes clear that without the Coulomb 
interaction U the pole €4 would be the lowest energy. According to the spectral weight a^a- 
in eq. f|T7|) . however, £4 requires a double occupancy of the lattice site. €4 to be the lowest 
energy therefore appears unphysical.On the other hand, for strong enough U it will not play 
any role. 

It is easy to generalize this alloy analogy to a disordered KLM. We have to take into 
consideration the non-magnetic sites €5 = as the fifth alloy constituent with the spectral 
weight 1 — X. The excitation spectrum is then: 

(19) 

£5 = -> a^^ = 1 - X. 

This zero-bandwidth alloy analogy is the basic framework for applying CPA to get the 
electronic selfenergy. 

5 

^AE) = E"-^i_(,^_s,(^))G„(E_s.(^))' (20) 

where 

G,{E) = ^y-^. (21) 
^ ' N^E-t^ ^ ^ 

k 




FIG. 2: Density of electron states for the disordered KLM at a)x = 0.95;b)x = 0.1, U — > oo, W = 
O.SeV, n = 0.06, S = 5/2, J = OAeV for ferromagnetic (S^) = S and paramagnetic (S^) = cases 

Therewith we can write the electron Green function: 

(E) = ^- , (22) 

where tj: is the Fourier transform of the hopping integral tij. The quasiparticle density 
of states (DOS) is derived essentially from imaginary part of the Green function: 

P^^E) = -^Y.^mG^K^E). (23) 

k 

Fig. [2] represents the quasiparticle density of states for large and small concentrations 
X of M atoms in full saturation (5^) = S and paramagnetic (S*^) = limits, respectively. 
In the case of strong Coulomb interaction(t/ oo), the DOS consists in general of three 
subbands. For (5*^) = S{T = 0) the spin up electron density is absent for energies around 
J{S + l)/2, while the spin down density is finite there. 

In the limit x = 0, our model reduces to the well-known Hubbard model. It has been 
shown j2l| that in infinite lattice dimension the CPA is an exact treatment of the alloy 
problem. But, the question, however, is: what is the optimal alloy analogy for the Hub- 
bard model? The conventional atomic-limit starting point, which we have applied, may be 



questionable for the pure Hubbard model. It is known (Ref. 22|]) that this alloy analogy 



does not allow for describing band-ferromagnetism. But, since we are interested only in the 
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ferromagnetism due to the interband exchange J, we beheve that the alloy analogy (Eq. 14- 
18) is appropriate for the correlated Kondo-lattice model although the Hubbard limit can 
be improved, for e.g. by the modified alloy analogy |23j. 



B. Spin Subsystem 

Let us now consider a structurally disordered system of Nm = xN spins which can be 
described by the isotropic Heisenberg Hamiltonian 

H = -Y. J^,Xfxf [StS^^ + StSj) . (24) 

We use two-time temperature Green functions for the investigation of spin excitations. 
The Green function within Tyablikov approximation satisfies the equation of motion 

E{{St\S,))^ = 25,,Xt^{S-t)+ ^2^^ 
2 Y.jjLi Jij {{Xf (5*1) - (S'f ) Xl^'^S^lS^))^ . 
After Fourier transformation to wave vectors, the above equation of motion can be written 

as: 

{E-x{S^)Eoiq)){{s^\S^,))^ = 

2x (S^) 6{q + q') + 2 (S^) ^^X.^^,+ (26) 
where we have used 

^X^=X^-x^5{q), 

X,= ^ELxt^e-^^^^^ (27) 
Eoiq)=2x{S^){JiO)-Jiq)) 



The approximation {S^) = {Sf) is very often applied to disordered systems and means the 
neglect of structural fluctuations of spin positions. 

Averaging over all possible realizations of atomic configurations, the equation for averaged 
Green function can be written in the following form: 



(E - Eom {{Si\S^,))^ = 2x {S^)S{q + q')+ ^^^^ 



The equation contains a higher-order averaged Green function AXG. One can write the 
equation of motion for this function, multiplying it by AX and performing configurational 



average which will include terms like, AXAXG. In order to solve these equations, the 
following decoupling of configurational averages is used; 



where 



= ^i<l- k')^X._^k^X^,_.. (30) 

The equation exploits translation symmetry. 

Thus, for the averaged Green function, we eventually get; 



c+ic-\\ - I a'\ 2x{s^)+p{m 
Pi'i;E) = ^E,^^i^s{q-k), (31) 

where we have introduced the structure factor for the random distribution of magnetic atoms 



Siq) = [{S^)YAX,^AX 



1. Virtual Crystal Approach 



This is the simplest approximation for the magnon Green function. If we neglect all 
scattering processes (P = 0, S = 0) in Eq. (13T1) we obtain the following expression for the 
magnon Green function: 



(32) 



2. Low Quadratic Approximation 



The next simple approach for the structure factor S{k)^ while including some scatter- 
ing processes (P 7^ 0, S 7^ 0), is obtained using a cumulant method [42|, |43|. Using the 
definition ( 1271) of variables 



(33) 



and the cumulant formula for finding the average of the product of X* variables [42, |43|. 



XfXf =5,,P2(x) + [Pi(x) 



(34) 
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where P2{x) = x(l — x),Pi{x) = x are second and first cumulants, respectively. Finally, 
using the following approximation for the structure factor S{k), 

S{k) ^ S'" = - x)x, (35) 

we get the expression for the magnon Green function in the Low Quadratic Approxima- 
tion(LQA) [Ref. Q]. 



pCk;E) = ^%^ (36) 

n(h- - {J{k-q)-J{k)){J{k-q)-J{q)) 
^^'^^ ^> - E-Eoik) 

III. RESULTS 

Let us discuss some of the results obtained for finite temperature using the self-consistent 
calculation according to the procedure as sketched in Fig. [TJ 

Fig. [3] represents the distance-dependence of the effective (MRKKY) exchange integrals 
for different values of the interband exchange parameter J of the Kondo-lattice model. For 
the concentrated case (x = 1) and in the low coupling regime J — the MRKKY inter- 
action agrees with the conventional RKKY [sol. In this situation the effective exchange 
constant has a long-range character with strong oscillations in the real space. However, 
with increasing J the MRKKY interaction looses the long range character and transforms 
into a fairly short-range interaction (like double exchange) , where only the first few effective 



exchange parameters 39|] turn out to be important. We realize the same behavior for diluted 
systems x < 1. For better comparison with the conventional RKKY, Fig. [3] shows only the 
oscillation part of the MRKKY interaction J'^^^ / J"^. The conventional RKKY interaction is 
a continuous function of distance i?, while the MRKKY approach can provide only discrete 
results. The magnetic neighbors of a given magnetic ion are considered as ordered in 'shells '. 
The larger the shell number L, the larger is the distance from the given magnetic ion. The 
L-th shell is built up by the L-th nearest neighbors. Each point in Fig. [3] corresponds to 
a certain shell demonstrating the distance-dependence of the effective exchange parameter 
Jij^ ■ Generally, there are two contributions to the effective exchange interaction, namely a 
Kondo-scattering of electrons J by magnetic impurities and in addition the scattering of the 
electrons by the random distribution of magnetic atoms with the concentration x. It is not 
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FIG. 3: Dependence of the MRKKY interaction (Eq. [TOj) on the distance R between magnetic 
atoms in a simple cubic lattice (ii = 1 is a first shell(nearest-neighbors atoms), ii = \/2 is a second 
shell(next nearest-neighbors atoms) and so on) at fixed parameters U — > oo,W = 0.96^,5 = 
5/2, X = 0.1, T = 0, (5^) = 5 of the model for three different coupling J = 0.2eV, J = OAeV, J = 
0.9eV and three different values of a small electron concentration n/x = l,n/x = 0.6, n/x = 0.2. 
In the inset is the electron DOS for the same parameters. The Fermi edge is indicated for the band 
occupation n/x = 0.6. 
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FIG. 4: Dependence of the nearest-neighbor effective exchange integral Ji and the next-nearest- 
neighboring one J2 on the exchange couphng J at C/ — > oo,W = 0.5eV,n = 0.06, S = 5/2, T = 0, 
for different values of concentration x. 

possible to separate these two contributions since both of them play an equvally important 
role. 




0-5 1.0 1.5 2.0 O.O Q.5 1 O 1.5 2.0 

J.eV J,eV 



FIG. 5: The same as in Fig. H but for S = 1/2. 

Figs, m and [5] represent the dependencies of the nearest-neighbor and the next-nearest- 
neighbor effective exchange interaction (Ji, J2) on the interband exchange coupling J, and 



13 




FIG. 6: Electron DOS calculated with U ^ oo,W = O.beV, J = l.2eV, 5 = 1/2, n = 0.6 a) T = 
for different values of the concentration x; h) x = 0.9 for different values of the temperature T. 

for different concentrations x and for two different spin values S = 5/2 and S = 1/2. In 
the strong coupling regime (J > 0.8) of the concentrated system {x = 1) only the short- 
range interaction (Ji) is important, the other interactions are small (J2, J3,...) [39]. In 
the diluted and strong coupling case, however, the next-neighbor interaction Ji becomes 
smaller (Fig. Hj) than for x = 1. In the low coupling regime (J < 0.1) the results of the 
modified RKKY coincide with those of the conventional RKKY ( Ji ~ J^) (Fig. Hf5|) which 
is nearly x-independent. The situation is much more complicated and non-monotonic in 
the intermediate coupling regime. Dilution may even lead to an enhanced value for the 
nearest-neighbor localized moment interaction Ji compared with that for x = 1 (Fig. Hf5l) . 
The self-consistent calculation of the effective exchange parameters Ji, J2, J3, .. yields only 
a slight temperature-dependence of order 0.01 meV or even less. 

The electronic quasiparticle-DOS for different concentrations x is shown in Fig. [61 There 
are two parts with different physical meanings. One part consists of correlated bands cen- 
tered at —JS/2 and J{S + l)/2 due to the exchange interaction with magnetic atoms. The 
spectral weights of the two subbands are strongly temperature-dependent (see Fig. Et^b)). In 
the case of T = the spectral weight of the upper f subband disappears, but however, 
is present in the J, spectrum. The second part around = represents a non-correlated 
band being connected to the non- magnetic atoms (Fig. [2|6|) . This part of the spectrum is 
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FIG. 7: Magnon DOS calculated self-consistently by use of the VGA (a) and the Low Quadratic 
Approximation (b) with the same parameters as those in Fig. [6] for different values of the concen- 
tration X. 

practically temperature independent. 

In Fig. [7] results for the magnon density of states are shown. Fig. [Tl^a) shows the magnon 
DOS in the VGA approximation fl32l) for the saturated ferromagnetic ground state. We see a 
typical consequence of dilution: a shift to the low energy side with decreasing concentration 
X of magnetic atoms. The same is valid for the magnon DOS in the low quadratic approx- 
imation (Eq. [5^ but there are some marked differences (Fig. [7](b)). Usually, in the case of 
weak magnon interaction (Fig. [7](a)), the magnon DOS near zero energy can be expressed 
as DOS{E) = DE^, where D is a stiffness parameter j^. But for the strong magnon in- 
teraction, due to the disorder, the magnon DOS near zero energy has higher contributions, 
like as DOS{E) = DE^ + KE^, where K is the disorder parameter. These contributions 
are strongly modifying the magnon stiffness. For example, the huge low-energy part for low 
concentrations x can be observed only for strong magnon interaction (inset of Fig. [7|(b)). 
But in the VGA we have the normal DOS{E) = DE^ behavior (inset of Fig. C^a)). We 
conclude the same as in ref. j^, that the disorder is strongly influencing the low energy 
magnon DOS. 

After the self-consistent calculation for finite temperatures, we obtain magnetization and 
the resulting value of the Gurie temperature (Fig. [H](a), Fig. Etb)). In Fig. [Hl^a), the magne- 
tization curves are plotted which are obtained within the VGA for the magnon subsystem 
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FIG. 8: Self-consistent calculation of the dependence of the magnetization (S^) on temperature by 
use of the VGA (a) and the Low Quadratic Approximation (b) with the same parameters as those 
in Fig. [6] for different values of the concentration x. 

(Eq. [32]) . We see that the Curie temperature decreases with increasing dilution 1 — x. The 
same holds for the electron polarization (inset of Fig. [8]^a)): (ri-f — n^)/n. 

The magnetization curves in Fig. Mjo) are determined within the low quadratic approach 
(Eq. [35)1361) for the same parameters as for VGA in Fig. M^a). Since this approach includes 
magnon scattering processes more realistically than the VGA, the Gurie temperature for 
the same concentration x of magnetic ions is lesser than that of VGA. For small x, the 
low quadratic approximation predicts first order ferromagnetic-paramagnetic transitions 3] 
(dotted lines in Fig. Et^b)). 

The explanation of this fact is based on predictions of the different approximations for 
the low energy part of the magnon DOS (see Fig. H^a) j7](b)). There was no report on first 
order transition in Ref. |35i]. In our opinion this is related to the fact that spin disorder was 
treated in mean-field like approach, which is evidently comparable to our VGA results(see 
Fig. [8](a)), and which does not exhibit first order transition. 

IV. DISCUSSION AND SUMMARY 

There are many real materials where disorder plays an important role for electronic as well 
as magnetic properties (binary substitutional alloys, diluted magnetic semiconductors, per- 
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ovskite manganese oxides, spin glass materials, transition metal dielectrics, etc.). Of course, 
these real systems are much more complicated than what the simple Kondo-lattice model 
predicts (complicated crystal lattice, multi-band structure, hybridization effects, spin-orbit 
coupling). However, we believe that the main microscopic mechanisms are well described in 
terms of the characteristic KLM features. The final goal is to make a quantitative descrip- 
tion of those materials combining the present analytic model investigations with realistic 
'ab initio' calculations of the band structure as it was done previously for concentrated 



local-moment systems 
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4l|. 



In this paper we discussed the influence of moment disorder on the electron and spin 
excitation spectra of the random KLM. Starting from an alloy analogy based on the exactly 
known zero-bandwidth limit of the KLM we applied a CPA procedure to find out the reaction 
of the electronic spectrum on the random mixture of magnetic and nonmagnetic atoms 18|. 
The Hubbard term in the model Hamiltonian ([1]) helps to prevent states, which belong to 
double occupancies of lattice sites, to be ground states. The analytical expression for the 
electronic selfenergy has been used then to get the effective exchange integrals of the modified 
RKKY theory. The latter results from a mapping of the interband exchange onto an effective 
random Heisenberg model (Fig. [1]) which was subsequently treated in the spirit of the well- 
known Tyablikov approximation. The disorder in the localized spin system turned out to be 
the most involved part of our study (Fig. [1]). It was incorporated via the equation of motion 
method and the technique of configurational averaging. In order to decouple the higher-order 
averaged Green functions we used the approximation of independent fluctuations ( l29i) . The 
expression (1311) for the averaged magnon Green function is generalized by using the structure 
factor of disordered distribution of magnetic atoms over a crystal lattice. Here we also used 
an approximation fl35|) . identical to the low quadratic approach 16j. 

There is no direct interaction between the localized moments. Therefore, the collective 
order is caused by the indirect exchange interaction mediated by the itinerant band elec- 
trons. Consequently, the indirect momentum coupling strongly depends on electronic model 
parameters such as exchange coupling J and band occupation n (see Fig. [3]). A further 
important parameter is of course the concentration of magnetic atoms x. The interactions 
found by modified RKKY resemble to those of the conventional RKKY only in the low 
coupling limit (J ~ 0). In the large coupling regime the long-range interaction transforms 
into short-range interaction (see Fig. Hll5ll . where only the nearest-neighbor interaction is 
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important(Ji 7^ 0, J2 ~ 0, J3 ~ 0,...). But in the intermediate couplings we found that 
( Fig. H|5|) the effective exchange interaction between the locahzed magnetic moments is 
strongly non-linear for small concentration x of magnetic atoms. This effect was also con- 



sidered in ref. 
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291], and also recently discussed in ref. 
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36 



semiconductors. In ref. 



25, 



371] for diluted magnetic 



381 ] the influence of disorder on the RKKY interaction of two 



magnetic impurities was accounted for by a phenomenological damping factor. In general, 
this damping factor has to be calculated within the full Kondo-lattice model. 

Another important finding is an influence of the disorder on the magnon excitations 
for small concentrations of magnetic atoms (Fig. [Tj). We found rather different results for 
the magnetic excitations in, respectively, VGA (Eq. [32], Fig. [7|(a)) and the low quadratic 



approach 



161 ] (Eq. [33 Fig. [71(b)). It is clear that VGA yields too simple expressions for 



the magnon Green's function (Eq. [36]and Eq. [3T|) . Recently some treatments were proposed 

in the mean-field (VGA) approximation in order 
33| taking into account only the dilution, while 



where the Kondo-lattice model is considerec 



to explain real material properties 25 



32 



our results show that for small concentrations the disorder effects are also very important. 
Such effects play a crucial role for the understanding and controlling key-properties such as 
the Gurie temperature. We see that disorder changes the magnon-DOS very drastically in 
particular for low energy excitations (Fig. [7|(b)) which are decisive for the resulting values 



of the Gurie temperature 



,l7|. 



We also note that there is another method to calcu 



ate the exchange interactions in 
works in the limit of infinite 



ferromagnetic metals and alloys 4J]. This method 
magnon wavelength. The MRKKY interaction does not have such a constraint. However, 
it is possible that for the disordered case the MRKKY theory (Eq. [TOj) has to be improved. 
The same leads true correct for the Gallen equation (Eq. [T2l) . 

Some factors have not been included in our self-consistent calculation, for example, the 
environmental cluster effects for electron excitations. We are aware that this effect can mod- 
ify the magnetic properties. These are currently being under investigation. Furthermore, as 
mentioned at the end of Sect. 2.1, the incorporation of the Hubbard type Goulomb interac- 
tion can be improved. This will be done in a forthcoming investigation which aims at the 
competition between band-magnetism due to the Hubbard interaction and local-moment 
magnetism due to the interband exchange J of the Kondo-lattice model. The atomic-limit 
alloy analogy is then of course not an appropriate starting point. 



18 



In this paper we have restricted our considerations to the ferromagnetic ( J > 0) KLM, the 



standard "Kondo physics" (J < 0) 46|, |47|] is therefore excluded from the very beginning. 



The main goal was to work out the interplay between disorder and magnetic stability in 
diluted local-moment systems. We are aware that parts of our theory have to be refined to 
reproduce the standard "Kondo Physics" (J < 0). 
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